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Abstract 

Physical strands or sheets that can be modelled as curves or 
surfaces embedded in three dimensions are ubiquitous in nature, 
and are of fundamental importance in mathematics, physics, bi- 
ology and engineering. Often the physical interpretation dictates 
that self-avoidance should be enforced in the continuum model, 
i.e. finite energy configurations should not self-intersect. Current 
continuum models with self-avoidance frequently employ pairwise 
repulsive potentials, which are of necessity singular. Moreover 
the potentials do not have an intrinsic length scale appropriate 
for modelling the finite thickness of the physical systems. Here 
we develop a framework for modelling self-avoiding strands and 
sheets which avoids singularities, and which provides a way to in- 
troduce a thickness length scale. In our approach pairwise inter- 
action potentials are replaced by many-body potentials involving 
three or more points, and the radii of certain associated circles 
or spheres. Self-interaction energies based on these many-body 
potentials can be used to describe the statistical mechanics of 
self-interacting strands and sheets of finite thickness. 
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1 Introduction 

The physical world consists of interacting matter. Many systems arising in 
science can be modelled effectively with a finite number of distinct constituent 
particles or point masses, qi G M. 3 (i = 1, . . . , N). A rich variety of natural 
phenomena - such as the existence of distinct phases (i.e. solid, liquid, and 
gas), and transitions between them - may be understood merely on the basis 
of pairwise interactions in such systems, with the phenomenological pair- 
wise interaction potential being induced by more fundamental interactions 
at the atomic level. Many-body interactions between triplets or quadruplets 
of particles are usually only treated as a higher-order correction. When the 
particles are all identical, it suffices to introduce a single potential energy 
function V(a) that is defined for all scalar arguments a > 0. Then the po- 
tential between distinct particles with labels i and j is given by V(|<& — qj\), 
where | ■ | denotes the usual Euclidean norm. We will be particularly con- 
cerned with systems where the potential has a qualitative form akin to that 
depicted in Figure [H i.e. strongly repulsive or infinite at short distances, and 
with a well at some finite distance a w . Provided that the potential is suf- 
ficiently repulsive in the sense that V(a) — > oo as a — > 0, all finite energy 
configurations have distinct locations for all particles, i.e. qi ^ qj for all dis- 
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tinct i, j = 1, ... ,N. Consequently, in any finite energy configuration, the 
minimum over all pairwise distances \qi — qj\ yields a distance of closest ap- 
proach for the configuration that is positive. And it can be anticipated that 
the location a w of the potential well provides a characteristic length scale for 
this distance of closest approach. 

We shall consider here the case in which the interacting particles are 
not unconstrained, but are instead restricted to lie on, or close to, a D- 
dimensional manifold embedded in IR 3 , such as a curve (D = 1) or a surface 
(D = 2). Such systems are widely studied in many branches of science; exam- 
ples include polymers [I] and random surfaces or membranes 0,0- The 
potential that gives rise to the internal stresses that maintain the system close 
to a continuous manifold is not our concern here - there are many such teth- 
ering potentials for curves and surfaces [2] [H] that are entirely satisfactory. 
Rather our focus is on the self-interactions that can arise when the curve or 
surface is sufficiently deformed in M 3 so that very different parts of the man- 
ifold come together to form, or are close to forming, a self-intersection. Such 
self-interactions are not captured by the standard local tethering potentials, 
and an additional non-local self-interaction potential must be introduced. 
Our key result is that in order to have a singularity-free description in the 
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continuum limit, the self-interaction potential must depend on D + 2 or more 
points. That is, for curves a non-singular self-interaction potential must de- 
pend on three or more points (and such non-singular interaction potentials 
exist), while for surfaces non-singular self-interaction potentials must depend 
on at least a four-point interaction. Notice that within this framework the 
standard case of a discrete number of unconstrained particles corresponds to 
D = 0, so that pairwise interactions suffice. 

The usual descriptions of a self-avoiding curve or surface employ pairwise 
potentials that are singular and which lack a length scale that can be identi- 
fied with the physical thickness of the system. For example, a single uniform 
strand of cooked spaghetti of length L and thickness h would typically be 
modelled by a curve C, which may be interpreted as the centerline of the 
spaghetto, together with an effective potential energy [3], jl] 

E 2 [G] = [ L U(q(s),q'(s), . . .) dC s + v ^ ^ 6(q(s) - q(a)) dQ s dG a (1) 
Jo Jo Jo 

where U is a specified function, q(s) is a parameterization of the curve, dC 
is an element of arclength, 6 is the Dirac delta function, and v is a model 
parameter. The first term arises from tethering effects at the microscopic 
level, and yields effective internal stresses in the strand. The second term, 
with its singular pairwise potential, models self-avoidance of the strand. But 
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this model does not capture effects of the physical thickness parameter h, 
which indeed does not even appear. The analysis and numerical treatment 
of such singular models is plagued by divergences that can only be handled 
using sophisticated mathematical techniques jjj, 

In this article we develop models of self-avoiding curves and surfaces that 
are nonsingular, and which include an explicit thickness length scale. In 
particular, we replace the singular delta-function potential above by a smooth 
potential V(a) of the qualitative form shown in Figure [T] that is dependent 
upon a single scalar variable a. The only difference between our models for 
discrete sets of points, for curves, and for surfaces, will be our choice for the 
independent variable a. Moreover, our continuum models retain the following 
two desirable features of the unconstrained discrete case described above: i) 
all finite-energy configurations of the manifold are non-self-intersecting, with 
a positive distance of closest approach, and ii) the location a w of the well in 
the potential V(a) provides a characteristic length scale for the distance of 
closest approach of the manifold. The difficulty in constructing such models 
is that the standard choice of taking the argument a of the potential to be 
the Euclidean pairwise distance cannot satisfy our two desiderata when the 
underlying system is continuous. The reasons are discussed more precisely in 
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Section|2]below, but the basic idea is simple. For the curve depicted in Figure 
I2fa), one wishes to penalize true points of closest approach between distinct 
parts of the curve, as, for example, between points 1 and 4, without penalizing 
adjacent points from the same part of the curve, as for example points 1 and 
2. The pairwise Euclidean distance simply cannot distinguish between these 
two cases; in other words, it cannot distinguish between proximity of points 
that is forced by continuity of the manifold in any configuration, and the 
real phenomenon of interest, namely proximity of points due to large scale 
geometrical deformation in some configurations. 

Our principal result is that the above two desiderata can be achieved 
simply by taking the argument a to be a quantity other than the Euclidean 
distance between two points. Specifically, as discussed in Section |2J for the 
case of curves a can be taken to be the radius of the circle defined by three 
points. When the argument of V(ct) depends upon more than simply two 
points, for example triplets or quadruplets of points, we shall refer to V(a) 
as a many-body or multi-point potential. The use of many-body potentials is 
an essential ingredient in the models of continuous systems that we propose, 
and they should not be viewed as a higher-order correction to two-body or 
pairwise potentials, as is the case in discrete models. Indeed our proposal 
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for continuous models is to replace pairwise self-interaction potentials, which 
must be singular, with many-body self-interaction potentials, which need not 
be singular. For example, for the spaghetto problem, we suggest an effective 
energy of the form 



E 3 [G}= U(q(s),q'(s),...)de s + / / V(r(s, a, r)) dG s dC, dC T 
Jo Jo Jo Jo 

(2) 



where the last term in the more standard energy (^Q) has been replaced by a 
three-body potential V(r(s, a, r)) with V of the form described in Figure ^ 
Here r(s, a, r) is the radius of the circle defined by the three points s, a and r 
as discussed in Section |21 The repulsive potential in (j2J) can be finite because 
in the limit of three points coalescing along the curve, for example points 1, 2 
and 3 in Figure|2^a), the radius r tends to the local radius of curvature, which 
is well-defined and positive for twice continuously differentiable curves. (A 
more detailed discussion of smoothness assumptions is given in Section 2.) 
However, whenever different parts of a curve come together to form a self- 
intersection, there are triplets of points, such as 1, 2, and 4 in FigureEfa), for 
which the corresponding circle radius r approaches zero. The self-avoidance 
of the uniform spaghetto of thickness h is modelled by a hard-core potential V 
which is infinitely large when its argument is less than h and zero otherwise. 
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Likewise, as described in Section a certain generalized four-body potential 
can be used to model self-avoiding surfaces of finite thickness. 

It should be stressed that the difficulties we address simply do not arise 
in a many-body system with a discrete index for the particles. On the other 
hand, analytic treatments of interacting systems with a very large number 
of particles are often facilitated by making a continuum approximation, in 
which the discretely indexed particles qi are replaced by a field q(x) that 
is dependent upon a continuously varying independent variable x. The ap- 
propriate phenomenological interaction potential for the continuum descrip- 
tion is to be derived from the microscopic ones. In the case of continuous 
phase transitions, this is a powerful procedure because the critical behavior 
is unaffected by the precise microscopic interactions 0, [101, and the crit- 
ical exponents can be derived using field-theoretic techniques [7j, jHj. But 
it is the passage to the continuum limit which implies that any two-body 
self-interaction potential must be singular. In order to avoid such singulari- 
ties we suggest that in the continuum limit the effective potentials modelling 
non-local self-interactions should be many-body ones. 
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2 Self-interactions of curves 

We consider various distances, other than the usual Euclidean pairwise one, 
that can be associated with points on a curve. Here and throughout a curve 
C will mean a function q(s) £ 1R 3 of a variable (arc-length) s £ [0, L]. We 
shall consider only sufficiently smooth curves, specifically those that are twice 
continuously differentiable. This stands in contrast to some models in field 
theory where polymers are sometimes represented by curves that are contin- 
uous but not smooth, for example piecewise linear. In point of fact there is 
an emerging body of literature JI] JI] which suggests that our ideas could 
usefully be applied to a slightly larger class of curves, namely those with only 
a Lipschitz continuous first derivative, but we do not pursue such questions 
here. A curve will be called simple if it has no self-intersections; that is, if 
q(s) = q(<j) only when s = a. 

In it was shown that certain ideal shapes of knots are related to vari- 
ous physical properties of knotted DNA. Intuitively these ideal configurations 
can be described as having the property that for a given knot-type and pre- 
scribed length they are as far as possible from self-intersection. The idea of 
a three-point distance based on the radius of the associated circle was intro- 
duced in [HIl as one wa y t° make the notion of distance from self-intersection 
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mathematically precise. In ^7] the same three-point circular distance was 
used as an ingredient in the numerical study of the optimal shapes of com- 
pact strings. The properties and relations between all possible circular and 
spherical distance functions defined on curves are discussed at length in [TH] . 

In the present article we argue that these generalized circular and spher- 
ical distances also provide natural means for defining singularity free self- 
interaction energies of curves through a potential function V(a) with the 
qualitative form of Figure [T] that takes a multi-point distance as argument. 
For our purposes, a self- interaction energy will mean a functional E[G] that 
is finite for any simple curve C, and which tends to infinity as C tends to a 
non-simple curve. We remark that within the specific context of knot the- 
ory the construction of simple geometric self-interaction energies for curves 
has already received much attention; the case of a pure inverse power of a 
circular three-point distance was proposed in [16], and surveys of alternative 
approaches can be found in several Chapters of [TTi] . 

2.1 Two-point distance for curves 

Given an arbitrary simple curve C, the most intuitive approach to the con- 
struction of a scalar argument a for a self-interaction energy V(a) is to take 
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the usual pairwise or two-point distance function 



P(s,a) = \q(s) - q{a)\. 



In particular, a candidate energy functional E[Q] would then be the double 
integral 



with for example V(a) = a~ m . The basic idea is that, for m > 2, the 
integral tends to infinity as q(s) tends to q{o) with s/c, thus meeting the 
infinite-energy condition associated with self-intersections. However, such an 
integral is always divergent due to nearest-neighbor effects since (3(s,a) = 
when s = a, so that the energy E[G] is infinite for any curve C. 

To cure the above divergence problem one may consider regularizing the 
integrand by subtracting something equally divergent as s — > a, or molli- 
fying the integrand using a multiplicative factor that tends to zero at an 
appropriate rate as s —>■ a. (See for example ^H] and references therein.) 
In essence, these procedures introduce a length scale to compensate for the 
fact that there is no inherent small- distance cutoff for the pairwise distance 
between nearest-neighbors along a curve. Renormalization group techniques 
may then be used to extract critical behavior that is independent of this 
artificial cutoff length scale jl], [HI- 




12 



2.2 Three-point distance for curves 

An alternative approach to defining the argument a of the self-interaction 
energy is based on triples of points. To begin, consider any three distinct 
points x, y and z on a simple curve C. When these points are not collinear 
they define a triangle with sides of lengths \x — y\, \x — z\ and \y — z\, 
perimeter ^{x, y, z) and area A(x, y, z). Each of these quantities vanishes in 
any limit in which all the points coalesce into one, so they do not individually 
yield an appropriate length scale for self-interaction. On the other hand, 
certain combinations provide quantities that remain positive in coalescent 
limits. One natural combination is 



which can be identified as the radius of the circumcircle, i.e. the unique circle 
passing through x, y and z. By convention, we take this radius to be infinite 
when the points are collinear. 

Various properties of the three-point circumradius function r(x, y, z) 
were studied in [T^j. For our purposes we merely note that the domain 
of the function r(x,y,z) can be extended by continuous limits to all triples 
of points on C, distinct or not. For example, if x = q(s), y = q(cr) and 
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z = q{r) are three distinct points on C, then it is straightforward to show 
that 



where p(x) is the standard local radius of curvature of C at x. (Because we 
consider only curves that are twice continuously differentiable this coalescent 
limit exists at each point; see ^3] for further details when the underlying 
curve is not smooth.) From its geometrical interpretation we may also deduce 
that r(x, y, z) and its limits are invariant under translations and rotations of 
a curve. Moreover, whenever different parts of a curve come close to forming 
a self-intersection, there are points x and y for which the limits r(x, y, y) 
and r(y,x,x) are both equal to half of the distance of closest approach. 

A general class of self-interaction energies involving the circumradius 
function can now be defined. In particular, one may consider the energy 



r(s,a,r) = r{q{s),q(a),q(r)) 

and V is of the form introduced in Figure or more simply V(a) = a~ m 
with m an appropriately large exponent. In contrast to the pairwise distance 




(3) 




where 
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function j3(s, a), the circumradius function r(s, a, r) does not suffer nearest- 
neighbor effects due to curve continuity. In particular, E[C] is well-defined 
for any simple curve C. 

The condition V(a) —>■ oo as a — > simultaneously provides control over 
a curve at both a local and global level. For example, at a point s on the 
curve, the limit r(s, s, s) defined in (J3J) is just the local radius of curvature 
at that location and the potential V(r(s, s, s)) is finite as long as this radius 
is non-zero (or equivalently, as long as the local curvature is finite). In this 
case, the three-body potential V plays the role of a local curvature energy 
that encourages curve smoothness. On the other hand, whenever different 
neighborhoods of a curve come together to form a self-intersection, there are 
points of closest-approach (s, a) for which the limits r(s, a, a) and r(a, s, s) 
tend to zero, leading to infinite values for V(r(s, a, a)) and V(r(a, s, s)). 
Thus, the three-body potential is also a global self-interaction energy that 
acts to discourage self-intersections. 

Characteristic length scales for curve self-interactions can be identified 
depending on the functional form of V. For example, potentials with the 
general form introduced in Figure ^ provide a natural scale for modelling the 
steric self-interactions of material filaments with non-zero thickness. A single 
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generic potential (repulsive at short r, attractive for intermediate r) based 
on a three-body argument r(s, a, r) suffices to obtain both a swollen phase 
at high temperatures and a dense phase at low temperatures, along with a 
phase transition between them in analogy with the fluid-solid transition for 
unconstrained particles with a pairwise potential. In the more standard con- 
tinuum approach |2] |I] a similar transition is obtained by introducing both 
attractive singular two-body (in order to encourage collapse) and repulsive 
singular three-body (in order to account for self-avoidance) (^-function types 
of potentials. 

We remark that the forces derived from a three-body potential have a 
geometrical interpretation analogous to those for a two-body potential. In 
the two-body case, the force on each particle is directed along the line that 
contains the two particles. In the three-body case, the force on each of the 
three distinct particles is directed along a radial line through the particle and 
the center of the circle that contains all three of the particles (cf. Figure EJ). 

2.3 Other distances for curves 

The three-point circumradius function leads to a notion of distance between 
pairs of points that contains geometrical information in addition to the stan- 
dard distance. For example, if x = q(s), y = qr(cr) and z = q(r) are three 
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distinct points on C, then one can associate a distance to x and y according 
to 

pt(x,y) = \imr(x,y,z) = j X V . 

T^cr 2\SmU X yi\ 

where 9 X y> is the angle between the vector x — y ^ and the tangent vector 
to C at y. The function pt(x, y) can be identified as the radius of the unique 
circle through x that is tangent to C at y. Notice that the function pt(x, y) 
will typically not be symmetric since the circle through y and tangent at 
x need not have the same radius as the circle through x and tangent at 
y. We will refer to pt(cc, y) as the point-tangent function. As a matter of 
convention, we consider it to be a particular three-point function, because it 
corresponds to a limit of the three-point circumradius function. 

The point-tangent function pt(x, y) shares two important properties with 
the circumradius function r(x,y,z). First, the coalescent limit cr — * s is just 
the local radius of curvature, namely 

lim pt(x,y) = p(x). 

Second, whenever different neighborhoods of a curve come together to form 
a self-intersection, there are points x and y for which pt(x, y) and pt(y, x) 
are both equal to half of the distance of closest approach. As a consequence 
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of these two properties, a general class of self-interaction energies involving 
the point-tangent function as argument can be defined as 



where pt(s,cr) = pt(q(s) , q(a)) and V is a general potential function as be- 
fore. The energy E[G] is well-defined for any simple curve C. 

Other distances for curves based on four points, or limits thereof, can 
also be considered. For example, to any four distinct points w, x, y and z 
on a simple curve C one can associate the four-point distance R(w,x,y, z) 
defined as the radius of the smallest sphere that contains all four points. 
Usually there will be such a unique sphere, but if all four points happen 
to be co-circular, there are many spheres passing through them with the 
smallest having the radius of the circle through the four points. Any limit 
of the spherical radius function as the four points coalesce into one is always 
positive because it is greater than or equal to the radius of curvature of the 
curve at the limit point. If the coalescence point has non-zero torsion, one 
obtains the radius of the osculating sphere in the limit, see for example [201 
p. 25]. Moreover, just as for the three-point functions r and pt, whenever 
different parts of a curve come together to form a self-intersection, there 
are points for which R is equal to half of the distance of closest approach. 
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(But now the possibility of end-point effects must be explicitly excluded, as 
for example if the curve is closed.) Similar results hold for the symmetric 
tangent-tangent distance function tt(x, y) defined to be the radius of the 
smallest sphere containing x and y, that is also tangent to C at both these 
points. 

3 Self-interactions of surfaces 

Self-avoiding surfaces have been a subject of much interest and study in di- 
verse disciplines ranging from mathematics to biology (see for example [5], 
jO], EI] Our objective is therefore to extend to the case of surfaces, 

our construction of multi-point distances that can lead to non-singular re- 
pulsive energies that preclude self-intersection. Our first conclusion is that 
for surfaces the simplest two-point (Euclidean), three-point (circular), and 
four-point (spherical) distances all share a similar problem: they cannot dis- 
tinguish between proximity due to continuity and proximity due to geometry. 
We will demonstrate, however, that the notion of point-tangent distance for 
a curve has a natural counterpart for surfaces that can make this distinction, 
and which provides a suitable argument a for a self-interaction potential. 
Throughout our developments a surface § will mean a (twice continuously 
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differentiable) function p(u) G M 3 of a variable u G A C M 2 , and a surface 
will be called simple if it has no self-intersections; that is, p(u) = p(v) only 
when u = v. 

3.1 n-point distances for surfaces (n = 2,3,4) 

Just as for the case of curves, an intuitive approach to the construction 
of a self-interaction energy for a simple surface S is to consider a repul- 
sive potential dependent upon the pairwise or two-point distance function 
P(u,v) = \p(u) — p(v)\ as argument. As before, while (3(u,v) has the desir- 
able feature that it tends to zero as different neighborhoods of S approach a 
self-intersection, it also tends to zero as u — > v by continuity, and thus leads 
to singular interaction potentials. 

In contrast to the case of curves, one may show that for surfaces both 
the three-point function r(u,v,w) and the four-point function R(t,u,v,w) 
introduced earlier, also lead to singular self-interaction potentials. In partic- 
ular, these functions tend to zero by continuity in the coalescent limit. This 
conclusion may be established as follows. Consider any fixed point a; of § and 
let £ e (cc) be the sphere of radius e centered at x. For each sufficiently small 
e > the intersection S £ (x) n S of the sphere and surface is a curve. The 
radius r of the circle through any three points on this curve satisfies r < e 
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(because the circle lies on the sphere of radius e), and similarly the radius R 
of the smallest sphere through any four points on this curve satisfies R < e. 
Thus, by considering the limit e — > 0, we can find sequences of three distinct 
points in the surface for which r — > 0, and sequences of four distinct points 
for which R — > 0. 

3.2 Tangent-point distance for surfaces 

For surfaces, the argument of a non-singular self-interaction potential can 
be obtained by passing directly to tangent-point distances. In particular, to 
any two distinct points x and y of a simple surface S we may associate the 
distance 



where the vector n x is either of the two unit normals to the surface S at 
the point x. The function tp(x, y) can be identified as the radius of the 
unique sphere through y and tangent to S at x. When y happens to be in 
the tangent plane to § at x, the sphere itself degenerates into a plane and 
tp(x, y) becomes infinite. Notice that tp(x, y) need not be symmetric since 
the sphere through x and tangent at y need not have the same radius. We 
refer to tp(x, y) as the tangent-point function for surfaces. Moreover, we 




y-x 
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consider it to be a particular four-point function since the tangent plane to 
S at x may be constructed through a coalescent limit of three points. 

The tangent-point function tp(x, y) enjoys various properties analogous 
to those of its counterpart defined for curves. For example, consider any 
curve q(s) in S such that q(0) = x and let t x = <?'(0). Then 

h tp ^^ = K.^(o)| =p(x ' y 

where p(x, t x ) is the absolute value of the local normal radius of curvature 
to § at a; in the direction t x . Thus, in coalescent limits y — > x, tp(x, y) may 
assume limiting values between the maximum and minimum of p(x, t x ) over 
directions t x with x fixed. When different parts of § come together to form 
a self-intersection, there are points x and y for which tp(a;, y) and tp{y, x) 
are both equal to half of the distance of closest approach of the surface to 
itself. 

While tp(x,y) need not be continuous when y = x, one may unambigu- 
ously consider self-interaction energies for surfaces of the form 

E[S] = J J V(tp(u,v))dS u dS v 

where tp(u,v) = tp(p(it), p(v )) and V is a general potential function. As 
in the case of a spaghetto, the choice of a potential of the form introduced 
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in Figure [T] could be used to model a surface of finite thickness. Similar 
ideas could be used for the description of triangulated or discretized random 
surfaces. In particular, one may discretize the energy introduced above to 
obtain a self-interaction energy for the triangulated surface. Here three ver- 
tices of each triangle could be used to define a tangent plane, which could 
then be used to evaluate the tangent-point function tp. 

We remark that Mobius-invariant energies for self-avoiding surfaces in- 
volving pairs of spheres, each tangent to the surface, are discussed in |2~o] . 
However, energies based on the radius of a single such sphere are not dis- 
cussed. 

4 Summary 

Curves and surfaces with an associated finite thickness can be used as contin- 
uum models of strands (such as a spaghetto) and sheets (such as this page). 
We find that non-singular self-interaction energies for such self-avoiding con- 
tinuous systems cannot depend on the usual pairwise distance. However, by 
using certain many-body potentials, we are able to construct self-interaction 
energies which are non-singular and which include a mesoscopic length scale 
for the physical thickness. Each many-body potential is a function of the 
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radius of a suitably chosen circle or sphere. Our energies provide simple con- 
tinuum models of strand and sheet systems that can be used to study the 
equilibrium and non-equilibrium statistical mechanics of distinct phases at 
different temperatures, along with phase transitions. 

We conclude by highlighting an application of our work in the context 
of the protein problem. Small globular proteins are linear chains of amino 
acids which, under physiological conditions, fold rapidly and reproducibly in 
a cooperative manner into their native state conformations These con- 
formations are somewhat compact structures corresponding to the minima 
of an effective energy. Furthermore, for proteins, form determines function 
and, yet, the total number of distinct native state folds is believed to be 
only of the order of a few thousand j^Zj and are made up of helices, hair- 
pins or sheets. An important issue in the protein field is to elucidate the 
bare essentials that determine the novel phase adopted by biopolymers such 
as proteins. Unfortunately, polymer science, which is a mature and techno- 
logically important field, does not provide an answer to this question. The 
standard, simple model of a polymer chain is one of tethered hard spheres - 
in the continuum limit, self-avoidance in this model is captured by a singular 
pairwise potential as in the second term on the right hand side of equation 
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(JTJ). Such a model with an additional effective attraction arising from the hy- 
drophobicity or mutual aversion of certain amino acid residues to the solvent 
(water), fails to account for the protein native structure phase on several 
counts. First, a generic compact polymer phase has many conformations 
which neither provide for specificity nor are kinetically readily accessible. 
Proteins have a limited number of folds to choose from for their native state 
structure and the energy landscape is vastly simpler. Second, the structures 
in the polymer phase are not especially sensitive to perturbations and are 
thus not flexible and versatile as protein native state structures are in order 
to accommodate the dizzying array of functions that proteins perform. 

Recent work [2E], built on the ideas presented here, has pointed out a 
crucial missing ingredient in the simple model of a polymer chain of teth- 
ered hard spheres. Strikingly, the simple physical idea of a chain viewed, 
instead, as a tube leads to several dramatic consequences. We are condi- 
tioned to think of objects as spheres and the effective interactions between 
them as being pair- wise in nature. This bias in our thinking arises from our 
everyday experience with unconstrained objects. However, when one deals 
with objects tethered together in a chain, some of the old notions need to 
be discarded. For example, a pairwise interaction only provides information 
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regarding the mutual distance between two interacting particles but does not 
have any contextual information regarding how far apart the two particles 
are along the chain. Our work here shows that in order to capture the con- 
straints imposed by a tube geometry associated with a discrete chain, the 
conventional notion of pairwise interactions between particles has to be aug- 
mented by appropriately chosen three-body interactions to capture the steric 
constraints imposed by the tube. 

Such a tube-like description of a chain leads to many of the standard 
results of polymer physics when the tube thickness is small compared to 
other length scales in the problem [21] • However, when the tube size becomes 
comparable to the range of the effective attractive interactions resulting from 
the hydrophobicity, the novel phase populated by biopolymers results |2*%] . 
This self-tuning of length scales occurs naturally for proteins because the 
steric effects (2H1 ED] associated with the backbone and side chain atoms 
of the amino acids, on the one hand, lead to a tube like description and 
control the tube thickness and the same side chains, on the other hand, have 
an effective attractive interaction which is at an Angstrom scale due to the 
screening effects of the water. 

It has been demonstrated [2H] that this novel phase has several charac- 
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teristics of the phase populated by protein native structures including the 
ability to expel water from the hydrophobic core, a vast simplification in the 
energy landscape with relatively few putative native state structures, a pre- 
diction that helices, zig-zag hairpins and sheets of the correct geometry are 
the building blocks of protein structures, a simple explanation of the cooper- 
ative nature of the folding transition of globular proteins and an explanation 
of why protein structures are flexible and versatile. Thus the idea of a tube 
and the use of appropriate many-body potentials introduced here not only 
lead to a better description of a polymer but allows one to bridge the gap 
between polymer science and protein science and provides, for the first time, 
a framework based on geometry for an understanding of the common char- 
acter of all globular proteins |2E] • We look forward to similar applications of 
our new ideas on sheets of non-zero thickness. 
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Figure 1: Example potential function with hard-core and potential well 
parameters a c and a w . 
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Figure 2: Interpretation of multi-point distances for a curve and a surface. 

(a) , Three-point distance for curve. Given any three distinct points, r is the 
radius of the unique circle that contains the points. When the points are 
from the same neighborhood on a curve, such as points 1, 2 and 3, r is close 
to the local radius of curvature. When points, such as 1, 2 and 4, are taken 
from two different neighborhoods of the curve that are close to intersection, r 
approximates (half of) the distance of closest approach of the curve to itself. 

(b) , Tangent-point distance for a surface. Given two distinct points 1 and 2 
on a surface, tp is the radius of the unique sphere that contains both points 
and is tangent to the surface at point 1. When the points are neighbors on 
the surface, tp approximates the absolute value of the local normal radius of 
curvature in the direction defined by the two points (not illustrated). When 
points are taken from different neighborhoods that are close to intersection, 
tp approximates (half of) the distance of closest approach. 
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Figure 3: Graphical interpretation of forces derived from potentials depend- 
ing on two- and three-point distance functions, (a), two particles interacting 
via a potential depending on two-point or pairwise distance. Holding particle 
2 fixed and moving particle 1 on the circumference of the circle leaves the 
energy unchanged. This implies that the force on each particle is along the 
joining line, and that the resultant of the forces is zero due to the translational 
invariance of the energy, (b), three particles interacting via a potential de- 
pending on three-point distance. Here any combination of the particles may 
be moved along the circumference of the circle without changing their inter- 
action energy. This implies that for any distinct triplet the force on each 
particle is along a radial line from the center of the circle to the particle, 
and that the resultant of the forces is again zero due to the translational 
invariance of the energy. 
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